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Abstract. We investigate the thermodynamics of artificial square spin ice systems 
assuming only dipolar interactions among the islands that compose the array. The 
emphasis is given on the effects of the temperature on the elementary excitations 
(magnetic monopoles and their strings) . By using Monte Carlo techniques we calculate 
the specific heat, the density of poles and their average separation as functions of 
temperature. The specific heat and average separation between monopoles with 
opposite charges exhibit a sharp peak and a local maximum, respectively, at the 
same temperature, Tp « 7.2D/kB (here, D is the strength of the dipolar interaction 
and ks is the Boltzmann constant). As the lattice size is increased, the amplitude 
of these features also increases but very slowly. Really, the specific heat and the 
maximum in the average separation dmax between oppositely charged monopoles 
increase logarithmically with the system size, indicating that completely isolated 
charges could be found only at the thermodynamic limit. In general, the results 
obtained here suggest that, for temperatures T > Tp, these systems may exhibit a 
phase with separated monopoles, although the quantity dmax should not be larger 
than a few lattice spacings for viable artificial materials. 



PACS numbers: 75.75.-c, 75.60. Ch, 75.60.Jk 
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1. Introduction 

New methods for exploring geometric frustrations in magnetic systems have been 
developed recently. Such methods consist in creating arrays of nanomagnets designed 
to resemble the disordered magnetic state known as spin ice. They are essentially 
composed of lithographically defined two-dimensional {2d) ferromagnetic nanostructures 
(elongated permalloy nanoparticles) with single-domain elements organized in diverse 
types of geometries (square lattice [U, hexagonal, brickwork [|], kagome ^ etc). 
Since their geometries are determined lithographically, lattice symmetry and topology 
can be directly controlled, allowing experimental investigation of a vast set of important 
theoretical models of statistical physics 0. These artificial magnetic compounds have 
the potential of increasing our understanding of disordered matter and may also lead 
to new technologies. Therefore, artificial spin ices are object of intense theoretical and 
experimental investigations H, |, |, |, |, 0, |, |, 0, [U], |13|, |14 . 



The trouble is that, in artificial spin ice patterns, the magnetization is unaffected 
by thermal fluctuations because the magnetic islands contain a large number of spins. 
Despite the fact that the moment configuration is athermal, these artificial materials can 



be described through an effective thermodynamics formalism |T6| ; in addition, some 
works have introduced a predictive notion of effective temperature 0, |T6[. For instance, 
an external drive, in the form of an agitating magnetic field behaves as a thermal bath 
and controls the temperature 0, [T^. Alternatively, this problem was addressed very 



recently by using a material with an ordering temperature near room temperature [T7 
such experimental work on a square lattice in an external magnetic field confirms a 
dynamical "pre- melting" of the artificial spin ice structure at a temperature well below 
the intrinsic ordering temperature of the island material, creating a spin ice array that 
has real thermal dynamics of the artificial spins over an extended temperature range 
||17|| . These findings and also other future possibilities make evident that a more detailed 
analysis of the effects of thermal fluctuations on a lower dimensional spin ice material 
should be of large interest for a better understanding of these frustrated systems. In 
particular, it would also be important to know the roles of elementary excitations in the 
thermodynamic properties of artificial magnetic ices. 

The main aim of this work is exactly this investigation. We are interested in the 
temperature effects on the excitations ("magnetic monopole defects" and their strings). 
Actually, since the prediction of monopoles in the usual three-dimensional {3d) spin 
ice materials |T8| and their experimental detection |19|, ^ |22|, the search for 



these objects in artificial compounds has become an important issue 11, 10, ^. The 
possible existence of these excitations in artificial and controllable systems is of great 
interest because they could be studied at room temperature and, more important, they 
could be directly observed with modern experimental techniques. Curiously, in the case 
of artificial systems, while the square lattice was the first to be produced |jl|, the direct 
observation of magnetic monopole defects and their motion was firstly accomplished in a 
kagome geometry Still, in this kagome lattice, a direct, real-space observation of the 
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Figure 1. (Color online) Specific heat as a function of temperature. It exhibits 
a sharp peak, at a temperature Tp ~ 7.2D/kB, which the amplitude increases very 
slowly with the system size L. Inset: the specific heat peak diverges logarithmically 
with the system size L. 



interplay of Dirac strings and monopoles was reported by Mengotti et.al For a square 
lattice, the direct observation of such excitations came only afterward because there was 
a primary experimental problem: until last year, none of such systems had achieved its 
ground state through thermodynamic equilibrium [|l^. Despite predictions [|, |, |^, 



the studies till recently do not have shown a long-range ordered configuration, perhaps 
because the researchers have used only non-thermal methods to randomize the array. 
This problem was experimentally solved by Morgan et. al. [|lOl- These authors have 
reported that by allowing the magnetic islands to interact as they are gradually formed at 
room temperature, the artificial square spin ice can be effectively thermalized, allowing 
it to find its predicted ground state very closely; thus, they could also identify the 
small departures from the ground state as elementary excitations of the system, at 
frequencies that follow a Boltzmann law. Subsequently, Magnetic Force Microscopy 
(MFM) images of a large number of isolated excitations with their string shapes and 
corresponding moment flip maps were described in square lattices |]10|. Therefore, 
the experimental results considering magnetic artificial square ices obtained in Ref . [p!0[] 
(which demonstrates the thermal ground-state ordering and the elementary excitations) 
and Ref. (which achieves a thermodynamic melting transition by using a material 
with ordering temperature near room temperature) lead us to think that more progress 
on the development of such arrays may become available in the near future, establishing 
opportunities to experimentally elucidate their real thermodynamics. 

2. The model and outlook 

Here, we consider an arrangement of dipoles similar to that experimentally investigated 
in Ref. [|l|. In our approach, however, the magnetic moment ("spin") of the island 
is replaced by an Ising-like point dipole at its center. In this approach, the internal 
degrees of freedom of each island are not being considered, as well as higher order 
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interactions. We expect that this simphfication does not change significantly the main 



physical properties of the system. As shown in Ref. ||2^, if the lattice spacing is about 
two times larger than the island's longest axis, the effect of higher order interactions is 
negligible. For smaller lattice spacings the effect of higher order interactions is to give 
more stability for the lowest energy states. In this way one may expect that as the island 
size increases, approaching the lattice spacing, the ground-state should be more robust 
and the appearance of excitations would cost more energy. While the consideration 
of the internal degrees of freedom would reduce the energy scale, the consideration of 
higher order interactions would increase it, but none of them are expected to change the 
physical picture discussed here. Thus, in our approach, at each site (xj, yi) of the square 
lattice two spin variables are defined: Sx(i) with components = ±1, 5*^ = 0, 5*^ = 
located at rx = {xi + 1/2, yi), and Sy(i) with components Sx = 0, Sy = itl,^^ = at 
fy = {xi,yi + 1/2). Therefore, in a lattice of volume = t^a^ (a is the lattice spacing) 
one gets 2 x P spins. Representing the spins of the islands by Si, which can assume 
either Sx(i) or 5*^(4), then the artificial spin ice is described by the following Hamiltonian 

Si ■ Sj 3(5*4 ■ ' '^ij 



Hsi = Da'Y. 



ij ij 



(1) 



where D = fiQfi'^ /Ana^ is the coupling constant of the dipolar interaction. We perform 
standard Monte Carlo techniques to obtain thermodynamic averages of the system 
defined by Hamiltonian (|I|). Periodic boundary conditions were implemented by means 
of the Ewald Summation EBI, used here to avoid spurious results brought about by 



the use of a cut-off radius|^^. Our Monte Carlo procedure comprises a combination of 
single spin fiips and loop moves where all spins contained in a closed random loop 
are fiipped according to the Metropolis prescription. In our scheme one Monte Carlo 
step (MCS) consists of 2 x single spin fiips and 0.7 x P worm moves. Usually, 10^ 
MCS were shown to be sufficient to reach equilibrium configurations and we have used 
10^ configurations to get thermodynamic averages. 

Before presenting the Monte Carlo calculations, it would be interesting to remark on 



some previous results [|, |TT|, ^ and some expectations for these arrays. The ground state 
configuration of the system in a square lattice is twofold degenerate. If one considers 
the vorticity in each plaquette, assigning a variable a = +1 and —1 to clockwise and 
anticlockwise vorticities respectively, the ground state looks like a checkerboard, with 



an antiferromagnetic arrangement of the a variable 10]. Of course, the ground state 
clearly obeys the ice rule (two spins point inward and two point outward in each vertex) , 
but with configurations of topology 1 (in 2d, there are two topologies that obey the ice 
rule. However, they are not degenerate and topology 2 is more energetic than topology 
1; see Refs.|Q Q for more details). The most elementary excitation is related to the 
inversion of a single spin (dipole) to generate a localized pair of defects. This is the 
3 — in, 1 — out state in a particular vertex and the 3 — out, 1 — in state in its adjacent 
vertex. In principle, these defects could be separated without further violation of the 
ice rule. Indeed, in our previous papers P, |^, we have numerically shown that these 
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Figure 3. (Color online) The average separation between charges exhibits a 
maximum around the same temperature Tp in which the specific heat has a sharp 
feature. The inset shows, in more details, the region around the maximum. 



defects behave as a monopoles pair since their interaction follows a d = 3 Coulomb law 
q/R, where q measures the strength of the interaction and R is the distance between 
the poles. However, we have also pointed out that an isolated monopole should be 
hard to see as effective low-energy degrees of freedom in the 2d square spin ice because 
the background antiferromagnetic order in the ground state confines them p|, since 
the ice rule is not degenerate in two dimensions. Actually, in 2d, there are additional 
excitations not present in the usual 3d spin ice [|1^], namely, energetic one-dimensional 
strings of dipoles (resultant spins at each vertex along a line of adjacent vertices) that 
terminate in monopoles with opposite charges. Such string excitations could be seen as 
lines which pass by adjacent vertices that obey the ice rule but sustaining topology 2 
(instead of topology 1) and hence they cost an energy equal to b times their length X, 
where b is the string tension. When the temperature T of the system is near absolute 
zero, the shortest path length connecting the monopoles gives the potential energy. The 
most general expression for the total cost of a pair of monopoles separated by a distance 
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R is the sum of the usual Coulombic term roughly equal to q/R, and a term roughly 
equal to bX resulting from the string joining the monopoles (there is, of course, also a 
constant term associated with the creation energy of a pair). Note that there is not a 
unique identification of a given path connecting the ends (monopoles) of the excitation. 
It is explicitly considered in the fact that the energy is proportional to X, which can 
assume different values for a given R. For a sufficiently long string, the string energy 
is completely dominant; for a short string the Coulomb interaction may have some 
importance if the size of the end-point monopoles is even smaller (as always occur for 
these systems). With the above features, these excitations are, to some extent, more 
similar to Nambu monopoles than Dirac monopoles. Really, as Nambu suggested, 
for a modified Dirac monopole theory, the string connecting monopoles has energy and 
is oriented, having a sense of polarization . 



In the artificial square ices, the ordering causes an anisotropy in the system making 
the monopoles interaction highly dependent on the direction in which the monopoles are 
separated in the crystal plane 0. This anisotropy is manifested in both the Coulomb 
and linear terms of the potential in such a way that we explicitly write 

V{R) =q{(^)/R + b{4>)X + c (2) 



where is the angle that the line joining the monopole defects makes with the x- 
axis of the array. Numerically, for instance, g(0) ~ —3.88Da, b{0) ~ 9.8D/a while 
g(7r/3) ~ —4:.lDa, &(vr/3) ~ 10. ID/a. The constant c ~ 23D, associated with the 
pair creation energy |0 {Ec ~ 29D) is independent of 0. Similar results can be found 
in the experimental work for the square lattice. Indeed, in Ref. [10 , the authors have 
found that, at a temperature T, these excitations arise in the system according to the 
Boltzmann law ~ exp{-/3V{R)) with b ^ WD /a, V{a) = ^ 30D and (3 = l/fc^T, 
where ks is the Boltzmann constant. They have also classified the elementary excitations 
by the number of fiipped spins (given by n) and a mnemonic character for shape. The 
three most observed defects are represented by 1 (a single pair with charges separated 
by only one lattice spacing) followed by 2L (a pair with n = 2 with the shape of L) 
and 40 (an isolated string loop with no charges and having n = 4 fiipped spins) ||T0[| . 
Curiously, the second excited state should be 40 since its energy is smaller than the 
energy of 2L defect. 

In principle, for the thermodynamics of these systems, the following argument 
should be valid: at low temperatures, there is insufficient thermal energy to create 
long strings (with length X larger than one lattice spacing) and so, the monopoles 
(with opposite charges) are bound together tightly in pairs. On the other hand, as 
the temperature is increased, the average separation between the constituents of a pair 
should also increase, which means that larger strings may become present in the system. 
Of course, there are several ways of connecting two monopoles by a string of length 
X. Therefore, considering states with X » R, we remember then that the number 
of configurations for the m-step self-avoiding random walk is = 5™, where 5 is a 
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Figure 4. (Color online) The maximum of the average separation dmax between 
opposite charges increases logarithmically with the system size L. 
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Figure 5. (Color online) Snapshot of a particular configuration of excitations for 
a temperature T = Q.OD/ks in a lattice with L = 10a. Red and black circles are 
positive and negative charges respectively. In general, for all temperatures below Tp, 
each monopole is clearly confined to its counterpart by a string (see the blue arrow 
indicating the direction of the string for the larger pair. Small pairs (i.e., monopoles 
bound together tightly in pairs) are indicated by a green arrow. 

constant and equal to 3 for a 2d square lattice. For the string with sufficient large X, 
N is well approximated by the random walk result and one obtains ~ 5^/". So the 
entropy of strings is proportional to X, i.e., the many possible ways of connecting two 
monopoles with a string give rise to a string configurational entropy proportional to X. 
Crudely speaking, then, the string free energy F = [b — {ln3)kBT/a]X will imply in 
an effective string tension [b — {ln3)kBT/a] which is positive in the low temperature 
region and the monopoles are completely confined. Above a certain temperature, it 
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Figure 6. (Color online) Snapshot of a particular configuration of excitations for a 
temperature T = 7 .&D /ks in a lattice with L = 10a. Red and black circles are positive 
and negative charges respectively. For a temperature above Tp, a small amount of 
monopoles does not have a string connecting them to their counterparts and, therefore, 
they seem to be isolated. There are also some pieces of strings (i.e., one-dimensional 
regions obeying topology 2, as indicated by blue paths) that do not connect monopoles. 
Small pairs are indicated by a green arrow. 

becomes negative, namely, the string looses its tension. The tension decreases like 
\b — {\n.?))kBT / a\ with increasing T, vanishing at some critical temperature /cbTc ~ 
6a/ln3. Using the average value for the string tension in Eq.(|^), i.e., h ~ 10Z)/a, we 
then estimate /cgTc ~ 9. ID. Of course, these theoretical arguments always overestimate 
the critical temperature. Although this picture leads to a rich physics for this system, 
predicting free magnetic monopoles and a phase transition, things may be a little 
more complicated. Really, additionally to the entropic effect discussed just above, 
there is another entropic contribution which manifests against monopole separation; 
the monopoles should become close together because it would provide more ways to 
arrange the surrounding dipoles in the lattice. Such effect introduces a 2d Coulombic 
interaction between the poles, which is proportional to T (i.e., K = Tln(i?/a)). If the 
temperature in which the string looses its tension is high enough, on the order of 9. ID 
as estimated, then, around this value of T, the confining potential Vs must be very 
strong, possibly preventing the freedom for the poles. With all these expectations, it 
would be important to investigate how the elementary excitations behave as a function 
of temperature. Our calculations is a first step in this direction. 
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Figure 7. (Color online) The density of string loops 40 (po) also exhibits a maximum 
around the temperature Tp ~ 7. 2D (green balls). This defect carries no charge and 
is the second excited state. Just for comparison, the density of pairs with opposite 
charges (ps) is also shown (red baUs). 

3. Results 

Now we present the results from Monte Carlo Simulations. The calculations shown here 
are for lattices with sizes 10,20,30,40,50,60 and 70 lattice spacings but in all figures 
we present only the results for lattice sizes 40, 60, 70. We start by presenting the results 
for the specific heat (see Fig.|I|). We notice that, for all lattice sizes studied, the specific 
heat exhibits a sharp feature at a temperature Tp approximately equal to 7.2D/kB- 
Indeed, the position of this peak does not seem to move as the lattice size L is varied. 
On the other hand, its amplitude Cmax increases much slowly as L increases. In the 
inset of Fig.|l|, we show how Cmax behaves with L. Therefore, with the obtained data we 
expect a logarithmic divergence of the specific heat in the thermodynamic limit. We also 
analyzed the pair density and the average separation between monopoles with opposite 
charg function of T. It is useful here to distinguish two types of monopoles: the 

less energetic ones in which the spins (in a vertex) are in the 3 — in, 1 — out or 3 — out, 
1 — in states (here referred to as unit-charged monopoles) and the most energetic ones 
in which the spins are in the 4 — in or A — out states (doubly charged monopoles). Figure 
H shows the density of pairs containing monopoles with unitary charge (ps) and also 
the density of pairs containing doubly charged monopoles {po, see the inset). They are 
calculated as the one-half of the thermodynamic average of the absolute value of the 
charge (±1) and (±2) respectively, summed over the lattice. For both cases, the density 
increases monotonously up to a maximum value achieved in the high-temperature limit. 

The size of the monopole pairs constitutes an internal degree of freedom, since the 
energy of a pair depends on the distance between the members of the pair. Here we 
would like to know the average distance rjv/ between two opposite poles as a function 
of temperature. Such a thermodynamic quantity may contain information about the 
possibility of monopoles separation and how they are organized into the system. For this 
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calculation we consider only defects with unitary charges. The grouping of monopoles 
into pairs is unique as long as the distances between them are smaller than the average 
distance between the monopoles r^/ = 1/ y/ps- As the size of the monopole pairs becomes 
larger than tm, one would simply have to redefine the monopole pairs. The average size 
tm of the monopole pairs is calculated by using the method of assignment problems; 
it deals with the question of how to assign n items (jobs, students) to n other items 
(machines, tasks) In our case, we would like to assign n positive charges to n 

negative charges for a given configuration in such a way that the sum of distances 
of all possible pairing be a minimum. The results are shown in Fig.^. The average 
separation has a local maximum at the same temperature Tp in which the specific 
heat exhibits a peak (~ 7.2D/kB). We notice that the amplitude of this maximum 
increases slowly as the system size increases. Indeed, like the specific heat peak, the 
maximum in the average separation dmax also increases logarithmically with the system 
size L {dmax oc InL, see Fig.^ and hence, one could expect that a certain quantity of 
monopoles may be almost isolated for very large arrays. Indeed, in our simulations for 
temperatures T >Tp considering lattices with L < 80a, we could observe some charges 
relatively distant from their respective counterparts (separated by distances of the order 
of 5a). For instance, we show in Fig. |^ a distribution of positive (red circles) and negative 
(black circles) monopoles in a small lattice with L = 10a observed in our simulations 
for a temperature T = Q.OD/ks (i.e., below Tp). Note that there are few excitations 
and all monopoles with opposite charges are coupled by a string, forming pairs. On the 
other hand. Fig. | shows the same system for a temperature above Tp (T = 7.6D/kB)- 
In this case, we see that a small quantity of monopoles are not connected by strings. 
In principle, they are free although some of them are not completely isolated (i.e., far 
away from other opposite poles). Furthermore, we also notice that some strings seem 
to be detached, not terminating in monopoles; there are few pieces of strings dispersed 
along the system (as said before, strings could be seen as lines which pass by adjacent 
vertices that obey the ice rule but sustaining topology 2 rather than topology 1). Of 
course, these figures exhibit only samples from a large number of data, but most of the 
data should be similar to the features of Fig.|^ for the regime of low temperatures and 
the features of Fig.|^ for the regime of high temperatures. Things must be clearer in the 
thermodynamic limit; in this case, some monopoles should become infinitely separated 
from their counterpart for temperatures T > 7.2D/kB- However, as the temperature is 
increased from zero, the monopole pair density grows simultaneously with an increase 
of the pair size (see also Fig. As the pairs become denser, there is less space to 
put in new pairs and hence the average pair size decreases for high temperatures. 
Really, we observe that, for T < Tp the average separation rM does not depends on the 
lattice size L, while for T > Tp, this quantity has a tiny dependence on L (at least in the 
range 7.2D /ks <T < 12D /kB). In this case, it is possible that monopoles may become 
completely isolated even for high temperatures {T > Tp) when L — )■ cxo. This picture for 
infinite systems corroborates the theoretical expectations for the existence of a phase 
with free monopoles |^ in large 2d artificial square ices, but the transition temperature 



Thermodynamics of elementary excitations in artificial magnetic square ice 



11 



( t t ( 
( ( ( ( 
( t ( t 


M ( t ) ( ( 
* t 1 ( ( M 

( * t ( ( t ( 


rTf 

1 t ( 




( t ( 
*rrr 
t t ( 


*rr 
t t 


( ( ) t 
*r(*rr 
( t ( ( 


t ( t ( 
( t ( ( 

1±H- 

( t ( ( 


* I t t t ( t 
t 1 f 1 t 1 t 
t t t t t ( ( 
( ( 1 ( ( t ( 

( ( t t t ( T 


( ( ( 
i t t 
( ( t 
( ( ( 
( ( ( 


( 1 1 1 1 ) 
1 1 1 1 1 1 
( ( 1 1 ( ( 
( 1 1 ) ( t 

M 1 t ( f 
( ( ( ) t t 
M 1 t ( t 


( ( ( 
( ( ) 
t t ( 


( * 
( ( 
( * 
( ( 
( ( 
1 ( 


MM 

t 1 t t 

t t \J 
( ( t ( 
t f M 


( 1 1 ( 
( ( ( t 
( t ( ( 

m 

fit:- 

( ( IT 
( t t t 


( * t ( ( t ( 


t ( t 
t t t 
( ( t 

i 


t * t t ( t 
( ( 1 1 ( ( 
t * t 1 ( t 

t ( t t ( r 

m 


( (Xrrr 
ICCCOX 

( t ( t ( n 
ICCCCCQ 


( M ( 
t ( M 
MM 

m 


1 ( 1 1 ( ( ( 

( ( 1 t T ( ( 

xcoxcc 

* ( » (X( ( 
t * t ( ( t ( 


( ( t 


* * T 1 1 ( 

"rrrrrr 
"rrrrrr 
( t ( ( ( t 




T 

4^ 

~r 
( 


t t t ( 
( ( ( t 

rrrt 


* t t ( t ) t 

:coxccc 

T ( ( T ( ) t 


( ( t 
t t t 
rrf 
rrf 
rrr 
rrr 
rrr 
rrr 
rrr 


( t t ) t ) 
t t t ) ( t 

:oxccc 

( t T t ( t 

"rrrrrr 
"rrrrrr 


( ( ( t t r 
t ( t 1 1 t 

~rrrrrr 
Trrrrr 
~rrrrrr 
Trrrrr 
~rrrrrr 


t ( 
( t 

~rr 
*rr 

-9r 


( 
( 

" 
*r 
~r 
"r 
~r 
"r 
~r 



Figure 8. (Color online) A typical configuration of string loops of the type 40 for 
a temperature below Tp (here, T — QD/ks)- At Tp, the number of 40 excitations 
proliferate in such way that a percolated cluster seems to be formed. The figure also 
shows the pairs of monopoles. 



(~ 7.2D/kB) should be little smaller than the estimated value ~ Q-lD/ks discussed 
earlier (remembering that the arguments of energy-entropy, in general, overestimate the 
correct quantity). 

We have also calculated the density of string loops 40, which is the defect with 
no charge but having the second lower energy (second excited state). Like the specific 
heat and the average separation, the density of defects 40 also displays a feature at 
Tp (see Fig. 0). Note that the string loops 40 almost do not appear in the system 
for temperatures smaller than Tp. Indeed, they surge suddenly at Tp and then, for 
temperatures above Tp, their number starts to decrease while the density of monopole 
pairs starts to increase more appreciable. Figures ( |) and ( P) show typical distributions 
of defects 40 in the system for temperatures below and above Tp, respectively. 



4. Discussion 



In summary, assuming the spin-spin interaction to be purely dipole-dipole, we notice 
that, at a temperature Tp, there is a maximum in the mean separation of opposite 
monopoles that increases logarithmically with the system size L {dmax oc InL). Hence, 
the distance between monopoles with opposite charges in the thermodynamic limit 
(T — )• oo) should diverge weakly, suggesting a possible unbinding of monopole pairs 
(T < Tp) into "free" monopoles (T > Tp). However, to the authors knowledge, for 
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Figure 9. (Color online) A typical configuration of string loops of the type 40 
for a temperature above Tp (here, T = SD/Zcs). The figure also shows the pairs of 
monopoles. 



a finite monopole density there is no diagnostic for (de) confinement based on a pair 
distribution function, for reasons analogous to tlie failure of the Wilson loop (which only 
knows perimeter laws in the presence of dynamical matter) to diagnose deconfinement 
in gauge theories. Indeed, from the three approaches that have been used to measure 
the static potential associated with the breaking of long flux tube between two quarks 
in QCD (i.e., correlation of Polyakov loops, variational ansatz and Wilson loops), string 
breaking has been seen only using the first two methods. On the other hand, the 
divergence found in tm could be understood in two different ways. It may be associated 
with either a vanishing string tension (which would lead to effectively free poles) or 
simply by the fact that in an order-disorder transition the correlation length (which is 
the only characteristic length of the system) diverges at the critical temperature. In 
this case, since the mean distance should be given in terms of the correlation length, 
then, it should also diverge. Of course, these two distinct ways to describe the system 
are closely related. We are faced thus with the question of the existence or not of a 
phase transition in this system. If there is a phase transition, other question arises: 
what is its nature? It is worthy to note at this point that although this system is 
closely related to the 16-vertex model, for which an exact solution is known, the range 
and symmetry of the interactions differ and thus we do not expect to observe the same 
critical behavior. Nevertheless, one point that deserves remark is the possible similarities 
between this system and the Ising model. In the two degenerated ground states, the 
a variables, related to the vorticity of each plaquette, can be seem as the spins of an 
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antiferromagnetic (AF) Ising model. In the AF Ising model, as the temperature raises, 
clusters of flipped spins are found in the system and at the critical temperature one can 
find percolated clusters of spins. If there is some similarities between these systems one 
may expect thus that the 40 excitations, which can be viewed as flipped a variables, form 
clusters at low temperature that percolates at the critical temperature, justifying thus 
the increasing number of these excitations at the transition temperature. This picture 
is corroborated by the logarithmic divergence of the specific heat. Unfortunately, our 
results are not conclusive about the possibility of a phase transition, and much more work 
has to be done in order to answer this question. To try to put some extra light on the 
topic, we have also done some calculations restricting the islands interaction to nearest 
neighbors converging in the same vertex, which would lead to a kind of generalized 
2c? Ising system with the same ground state. Nevertheless, we have obtained that the 
vertices with topology 3, in the 3-in/l-out and 3-out/l-in states, remain connected by 
strings (but now, there is no Coulomb interaction anymore). The interaction energy 
between two opposite vertices in topology 3 (type /// vertices) is given by bjX + cj, 
where bi = 26D/a and cj = SAD, much bigger than the usual results obtained for 
the long ranged dipolar interaction. Since the string tension remains, the arguments 
associated with the string configurational entropy should maintain valid and we have 
again the same problem as before (but with different energetics; for instance, the value of 
the temperature in which the quantities show a maximum changes to IGD/ks)- Indeed, 
the specific heat, the average separation between opposite type /// vertices etc, have 
the same behavior found for the system with long-range dipolar interaction (not shown 
here). 

From a practical point of view, the divergence in in the thermodynamic limit, 
and thus the phase of large separation among monopoles should not be expected in finite 
systems. Due to the slow logarithmic divergence, the extrapolation of our results to a 2d 
lattice containing the Avogadro's number {N^^^ = 10^^ = 10^ x 10*^) of islands will imply 
in dmax ~ 2.5a only. On the other hand, even with small values for dmax, some monopoles 
may become isolated for temperatures near 7. 2D (see Fig. P). The challenge of building 
arrays using new materials (with an ordering temperature near room temperature ) 
and/or with reduced island volume and moment (and possibly with larger L) should 
be then an important issue for technological applications. Indeed, it concerns with 
the excitations evolution in these artificial compounds. These developments may 
experimentally determine the possibility of monopole dynamics, their lifetimes and so 
on. For instance, based only on the average separation results, we speculate that, near 
the temperature Tp, the annihilation process of monopoles (without strings) should be 
more probable to occur in small arrays than in large arrays due to the fact that the 
mean separation between such opposite charges increases with the system size. 
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